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THE MAGNETIC COUPLING OF PLANETS AND SMALL BODIES WITH A 

PULSAR'S WIND 



F. Mottez^ and J. Heyvaerts^ 

Abstract. We investigate the electromagnetic interaction of a relativistic stellar wind with a planet or a 
smaller body in orbit around a pulsar. This may be relevant to objects such as PSR B1257+12 and PSR 
B1620-26 that are expected to hold a planetary system, or to pulsars with suspected asteroids or comets. 
Most models of pulsar winds predict that, albeit highly relativistic, they are slower than Alfven waves. In 
that case, a pair of stationary Alfven waves, called Alfven wings (AW), is expected to form on the sides of 
the planet. The wings expand far into the pulsar's wind and they could be strong sources of radio emissions. 
The Alfven wings would cause a significant drift over small bodies such as asteroids and comets. 
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1 Introduction 



Two pulsars with planets have been discovered, in 1992 and 1993 ( Wolszczan &: Frail 119921: iThorsett et al.lll99,Sh . 
In particular PSR B1257+12 hosts three planets at distances of the order of the astronomical unit. A pulsar is 
a neutron star with a fast spin (period P ~ 0.001 — 5 s) and surrounded by a magnetosphere that is a powerful 
source of radio and/or high energy emissions. (We use also the frequency 17* = 27r/P*). Precise timing of pulsars 
show that they spin down at a rate that is typically P ^ 10~^^ — 10~^^P. The dissipation of energy of PSR 
B1257+12 associated to the spin down is Erotation = -MjQM* = iir'^MiP/P^ = 2. x 10^'' W, where Mj is their 
momentum of inertia. It can be compared to the planet's gravitational energy Eo — GAI^,AIp/2a = 4. x 10'^^ 
J. If the flux dissipated by the spin down was captured through the planetary radius, the gravitational energy 
of the planet PSR 1257+12 "a" would be dissipated in only 8. x 10^ years. Even with an inefficient coupling 
mechanism, can we expect heating, plasma acceleration and radio emissions from the planet, or a modification 
of its orbit ? 



2 A planet in the wind 

The planets (with an orbital period Port > a few mn) are far beyond the light cylinder; they orbit in the 
pulsar's wind. This pulsar wind is an almost radial ex pansiori of ult r arelativistic (Loreiitz fac t or 7 ^ 10^ — 10^.) 

and underdense plasni a. Most pulsar wind models (Michel 1969t Contopoulos et al. 119991: ISpitkovskvl [20061 

Kirk et al. l l2009l:IPetrill20ld ) converge on the following facts: This wind is highly magnetized (B^ >> i^oPl), 
the magnetic field is quasi-azimutal {B >> Bpoioidai)- It is possibly sub-Alfvenic {vw ~ c, V^i ~ c, until 

lOO's of r^c, (models) vw < Va)- This last fact means that the planet is not preceded by a shock wave, but it 
is directly connected to the wind. 



3 An electromagnetic wake behind the planet: Alfven wings 



Then, making simple hypothesis, it is possible to develop a theory of the interaction of the p lanet with the wind. 
The hypothesis are (1) the plasi na is locally unif orm and incompres sible (OK at 1st order (IWright fc Schwartz 
1993)). It can be non-relativistic (|Neubaueij|l980t ) or relativistic flow ([Mottez fc Hevvaertsl2011bHHevvaerts et al 
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2012h . (2) Alfven waves (MHD waves carrying mainly magnetic fluctuations) propagate almost only along the 
magnetic field (3) the planet is a conducting body (maybe because of a ionosphere) 

Then, a potential drop is established by the plasma flow across the planet. The potential drop induces a 
current. This current is closed along the magnetic field lines, by two stationary Alfven waves called Alfven 
wings. This structure is a kind of electric circuit, sketched on see the left-hand side of Fig. [TJ Its conductivity 
along the planet is those of the planet. Its conductivity along the two ribbons of current that extend into 
space (the Alfven wings) is, in the highly relativistic case, approximately equal to the vacuum conductivity 



Y.A - l/Aioc--0.09 S. 



(3.1) 



Then, it is possible ( Neubauei 19801 : Mottez fc HevvaertsI 2011b ) to estimate the total current circulating in the 
Alfven wings 

Iaw = 4(£;o - E,)RpJ:a, (3.2) 

where Eq is the mean electric field associated to the unipolar inductor, Ei due to the planet's resistivity (that 
is low), Rp is the planetary radius, or more generally, the body's radius. For a P = 10ms PSR and an Earth- 
like planet at 0.2 UA, Iaw ^ 10^ A. It is larger for a standard pulsar (of period P = Is), Iaw ~ 10^^ A. 
For comparison, the pulsar's current i.e. the Goldreich- Julian current that powers the whole magnetosphere is 
Igj = 2. X 10^^ A. For standard pulsars, the two currents are comparable ! 



4 Radio emissions ? 

Of course, these strong and localized ribbons of electric current are probably unstable. Consequently, they 
may constitute powerful sources of radio emissions. A priori, there are two possible regions of emission (it is 
nothing more than a conjecture) (1) a radio source near the planet. Then, the source does not move in our 
frame of reference and the emissions might be similar to "usual" planetar y radio em issions, generally associated 
to instabilities w ith J||, as it happens with the couple Jupiter-Io (Que innec fc Z arka 1998; Hess ct al. 20^^ 
Hess et al. 20091) . It would be difficult to predict the angles of emission. (2) The radio source is along the 



Alfven wing, and it is convected with the pulsar's wind. Then, the source is moving fast V ^ c and the emission 
is beamed by the relativistic aberration (see the right-hand side of Fig. [l]). In that cas e, the radiation would 



be well collimated, and it might be observable, under lucky circumstances (iMottezll201ll) 



5 Alfven wings have a large influence on comets and asteroid's orbits 

In the Alfven wing electric circuit, the current density j fiowing along the planet, crossed with the ambient 
magnetic field, is the cause of a j x B force density. The force has a component that is tangential to the orbit. It 
gives a non Keplerian contribution t o the motion. This force and its effect on the orbit of a body in the pulsar 



wind can be estimated analytically (jMottez fc Hevvaertsll2011al) . Then, the evolution of the orbital elements 



can be computed numerically. For planets this force has a negligible contribution to the orbit. But it becomes 
very important when small bodies, like comets, asteroids and planetesimal are concerned. Figure [5] shows the 
time evolution of the semi-major axis of a 1 km sized rocky body in orbit around a standard pulsar (pulsar 
period P — Is). The forces acts differently if the orbital spin is parallel or anti-parallel to the neutron star' 
spin. When they are parallel, the semi-major axis and the eccentricity tend to increase, while they are reduced 
for anti-parallel spins. In this last case, we see that a body can be precipitated onto the neutron star in less 
than 10,000 years. 

This effect of the magnetic force on the objects orbiting a pulsar, especially planetesimals, may have impor- 
tant consequences in scenarios of planetary formation after the transformation of the star into a neutron star, 
from the debris of the supernova. 



6 Conclusions 

A planet orbiting around a pulsar would be immersed in an ultrarelativistic underdense plasma flow. 
It would behave as a unipolar inductor, with a significant potential drop along the planet. As for lo in Jupiter's 

magnetosphere, there would be two stationary Alfven waves (Alfven wings) attached to the planet. 
The AW are supported by strong electric currents, comparable to those of a pulsar. It would be a cause of 



planets and small bodies in a pulsar's wind 
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Name 


C/(V) 


/(A) 


Ejmax (W) 


/S.a/year (m) 


/S.e/year 


PSR 1257+12 a 


1.1 xlO^^ 


3.0 xlO^ 


2.5x10^^ 


0.02 





PSR 1257+12 b 


3,5 xl0i2 


9.4 xlO^ 


2.5x10^2 


1. x 10-3 


3 xlO-i'^ 


PSR 1257+12 c 


2,6 xl0^2 


7.0 xlO^ 


1.4x10^2 


1. x 10-3 


2.4x10-**^ 


PSR 1620-26 a 


6,0 xlO^^ 


1.5 xlO-^ 


7x10^" 


S.xlO-*^ 





PSR 10ms b 100 km 


2.4x10^ 


6.xl0'^ 


1.2x10^** 


8. x 10-3 


e.xio-^'' 


PSR 10ms b 1 km 


2.4x10^ 


6.x 10* 


1.2x10^2 


0.8 


6.x 10-^2 


PSR 1 s b 100 km 


2.4 xlO" 


6x10** 


1.2 xlO^" 


8x10* 


6.4 xlO-** 


PSR 1 s b 1 km 


2.4 xlO^ 


6.x 10*^ 


1.2x101*^ 


S.xlO^ 


6.x 10-*^ 



Table 1. Electric potential drop, total electric current associated to the Alfven wing. Electrical energy Ejmax dissipated 
in the Alfven wing. Variation per (terrestrial) year of the semi-major axis. Variation of the eccentricity, per year, 
Ae/year. 




Fig. 1. Left: The engine of the Alfven wing is a unipolar inductor. The unperturbed wind's magnetic field Bo and velocity 
vo are almost, but not exactly, perpendicular. The electric field Eo created by the unipolar inductor is perpendicular 
to these two vectors; it induces an electric current (of density j) along the body. This current then goes into the 
interplanetary medium, forming two structures, each of them made of an outwards and an inwards flow. The current 
density j flowing along the planet is the cause of a j x B force density. Right: general shape of the Alfven wings on a 
larger scale. The spiral is the shape of a magnetic field line. The arrows show the directions of emissions of the radio 
waves if they are emitted by a source that is propagated by the wind. 



strong radio emissions, with sources all along the AW, highly colimated through relativistic aberration. 
There would be a chance to detect these radio-emissions from Earth. 

The emission would be pulses as for ordinary pulsars, but highly dependent on the planet-star-observer angle, 
maybe one very brief sequence once every orbit. 
The Alfven wing exerts a force unpon the orbiting body. This force has no influence of the motion of a 

planet. On the time scale of millions of years, it can however affect the orbit of bodies with a diameter of 100 
kilometres around standard pulsars with a period P ^1 s and a magnetic field B ~ 10^ T. Kilometer-sized 
bodies experience drastic orbital changes on a time scale of 10* years. 
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Fig. 2. Time-evolution of the semi-major axis of a 1km sized asteroid under the influence of the magnetic thrust. The 
initial distance a is 0.16 UA. The different curves are traced for various values of the orbit eccentricity e. Left: The 
orbital spin is parallel to the neutron star spin axis. In that case, the semi-major axis increases. The time axis covers 
100,000 years. Left; The orbital spin is anti-parallel to the neutron star spin axis. The semi-major axis decreases. The 
time axis covers only 10,000 years. Right: 
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